No-soliton— soliton phase transition in a trapped one-dimensional Bose gas 
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Following the experimental observation of bright matter- wave solitons [L. Khaykovich et al, Sci- 
ence 296, 1290 (2002); K. E. Strecker et al, Nature (London) 417, 150 (2002)], we develop a 
semi-phenomenological theory for soliton thermodynamics and find the condensation temperature. 
Under a modified thermodynamic limit, the condensate occupation at the critical temperature un- 
dergoes a sudden jump to a nonzero value, indicating a discontinuous phase transition. Treating the 
condensation as a diffusion over a barrier shows that the condensation time is exponentially long as 
one approaches the thermodynamic limit, and the longest near the critical temperature. 
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Quantum-degenerate atomic gases are attracting at- 
tention as versatile physics laboratories with unprece- 
dented level of experimental control Nonlinear sci- 
ence in particular should benefit from the study of dark 
and bright 0, Q matter-wave solitons Q recently 
observed in one-dimensional (ID) Bose gases, whose 
mean-field dynamics is described by the ID nonlinear 
Schrodinger equation. Much recent theoretical work on 
ID matter-wave bright solitons is dedicated to explain- 
ing their formation, observed number, and dynamics 
and to proposing new ways of generating them 0. We 
shall consider instead the thermodynamics of this system, 
rarely of central importance in other nonlinear media. 

Microscopically, a bright soliton is a many-body well- 
localized bound state of an attractive gas. While at zero 
temperature this state contains all the atoms, at suffi- 
ciently high temperature it should disappear. We de- 
scribe the ID harmonically trapped gas using a model 
with exactly solvable thermodynamics and show that, 
in a modified thermodynamic limit, the change in the 
bound-state population happens through a discontinu- 
ous phase transition. Our main computational result is 
the critical temperature, with accompanying phase dia- 
gram. Finally, we turn to metastability issues and con- 
sider condensate formation (and evaporation, above the 
critical temperature) as a diffusion over a potential bar- 
rier. We find that the transition time increases roughly 
exponentially with a parameter characterizing the prox- 
imity to the thermodynamic limit, the duration peaking 
near the critical temperature. 

Ground state. The Hamiltonian for a ID system of N 
attractive (^-interacting bosons in a harmonic trap is 
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We assume the ID coupling constant to be negative, 
5id < 0; its relation to the ID scattering length <zid 
is <?id = — 2fi, 2 /(m oid). In the actual waveguide or cigar 
trap experiments the value of a\D is governed by the 3D 
scattering length a and size of the ground transverse state 
a± = ^/2h/muj±, namely am = — ( a l/2a) [l—C(a/a±)], 
where C = |C(1/2)| = 1.4603 . . . 0. Here u± is the fre- 
quency of the transverse confinement, m is the atomic 
mass, and £(x) is the Riemann zeta-function. 

In the absence of the trapping potential the ex- 
act ground state of the Hamiltonian assumes a 
known Bijl-Jastrow form 0, Hoi 111) ^(zi, . . . , 2jv) = 
^Yli<i<j<N e~\ Zj ~ Zi ^ aiD , where A is the normalization 
constant. Ground state energy is then £?Q Xact (iV) = 
-_!_ (mgjv/h 2 ) N(N 2 - 1) . The center of mass of a 
system in this state is delocalized, although the rela- 
tive motion of particles is tightly bound. If we now 
add a weak harmonic confinement Ut iap (z) = moj 2 z 2 /2 , 
(where the chemical potential of the free space soliton 
m(gN) 2 /h 2 3> hu> z ), the one-body density distribution 
becomes localized around the origin. In the limit of 
large N, mean-field theory predicts a one-body density 
distribution of the form p(z) = po / cosh 2 (z / £) , where 
po = (2^) _1 is the peak density and I = 2h 2 /mgN is the 
density profile width 0, [n| . 

Model Hamiltonian. We analyze the thermodynamics 
of our system using the following Hamiltonian: 



H = E {N 
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U(zi). (1) 



where E (N ) — —± (mgl^/h 2 ) N§ is the soliton en- 
ergy for Nq 3> 1, iVo = N — N' is the number of particles 
in the soliton, N is the fixed total number of particles, 
=o ^ S t ne number of noncondensed parti- 
cles, e s = hw z s is the one-body energy spectrum of the 
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ID harmonic oscillator, and uj z is the frequency charac- 
terizing the longitudinal confinement. The Fock vacuum 
for the bosonic annihilation and creation operators a s 
and aj, where [a s , ap] = 5 SiS i, is the harmonic oscillator 
ground state. Note that the ground state of the Hamil- 
tonian H corresponds to all N particles being in the soli- 
ton, so that the first excited state is N — 1 particles in 
the soliton, 1 particle in the lowest harmonic oscillator 
state s = 0. The energy spectrum thus has a gap, about 
equal to the chemical potential of the free space soliton, 
as can also be confirmed by Bogoliubov analysis [Toj ]. 

By far the most "phcnomenological" assumption we 
made is the conjecture that even at finite temperatures 
the condensate is not fragmented, but only depleted. 
This conjecture is confirmed in ab initio calculations for 
mesoscopic (N ~ 40) numbers of atoms in a ID box [l2j |. 
Next comes the assumption that the solitonic conden- 
sate is not deformed by the external harmonic confine- 
ment, leading to the requirement 01 ^ dEo/dNo ~ 
m(gNo) 2 /ti 2 . We further neglect interactions both be- 
tween the noncondensed particles and between noncon- 
densed particles and the soliton. The former can be justi- 
fied for the case of well-delocalized noncondensed cloud. 
To justify the latter notice that under the above require- 
ment of no deformation by the external confinement, the 
size of the soliton i becomes much less than the size of 
the trap ground state a z = y^2h/mu! z . One can show 
that strong but well-localized perturbation cannot affect 
the overall density of states. The assumption can also 
be justified via a Bogoliubov-like analysis of the prob- 
lem 0, 0] . Finally, in our Hamiltonian we neglect the 
center-of-mass motion of the soliton as it can be shown 
to be fully decoupled from the rest of the system. 

Finite size thermodynamics. Thermodynamics of our 
system is determined by the canonical partition function 



A" 



where 

Q'(N>, T) = 



Q'(N-N 0) T), 
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is the canonical partition function for TV' non-interacting 
bosons in a ID harmonic trap In Fig. da) we show 
the mean soliton occupation and its standard deviation 
due to thermodynamic fluctuations as a function of tem- 
perature for a typical experimental setup. In this calcu- 
lation we have replaced the sum over Nq by an integral 
and used a large temperature (k^T 3> huj z ) approxima- 
tion for the partition function of noncondensed particles, 
ft"* -yi 



Q'(N', T) 
where Li2(C) 



exp 

In(l-t) 



k B T 
huj z 



Li 2 (l) -Li 2 (e fc B T ) 



dt is the dilogarithm function. 



Scaling. The following scaling has proven useful 
for analyzing thermodynamical properties of our sys- 
tem. Let us introduce a dimensionless parameter e = 
(mg 2 D )/ (h 3 uj z ), a dimensionless rescaled temperature 
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FIG. 1: Solitonic condensate fraction vs. temperature for 
N — 1050 7 Li atoms of a scattering length 0315 = -1.46 nm, 
confined in a cigar-shape trap with frequencies u>± = 2% x 300 
Hz and ui z — 2n x 3 Hz. This set of parameters corresponds 
to a = 56 huj ± , e = 1.77 x 10" 4 , and 77 = 0.186. (a) Finite 
theory; (b) thermodynamic limit, e — > 0, with the value of 
the energy scale a and rescaled total number of particles rj 
kept constant. The critical temperature ksT c = 1.65 hui±_ 
corresponds to r c = .0293. 



r = knT/a, and a dimensionless rescaled number of par- 
ticles 77 = eN, where a = {h A u 2 ) / {mg^). We also in- 
troduce rescaled numbers of condensed (770) and noncon- 
densed {rj 1 ) particles in an analogous way. The parameter 
e can be shown to be proportional to the ratio between 
the binding energy of a ID dimer £(, = (mg 2 D ) / (4h 2 ) 
and the trap level spacing h\o z : e ~ (1/4) (eb/huj z ). The 
rescaled partition function is 



Q(r), t) 

feijl, t\t] ) 



sl dr] e 



-feh, r I Vo)/t 



±{ A% 3+t 2 [Li 2 (l)-Li 2 (e-^)]} 



Notice that the positions of the extrema of the Landau 
free energy / £ (?y, r 1 770) as a function of 770 does not de- 
pend on either the energy scale a or the parameter e. 
Here we attempt to identify the rescaled number of con- 
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Rescaled number of condensed atoms, T) = £ N 
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FIG. 2: Free energy as a function of the order parameter r? Rescaled temperature, x = k B T/c 

for the set of parameters corresponding to Fig. 1. 



densed particles t/q with the order parameter of a possible 
phase transition. The distribution of the order parameter 
is given by the Boltzmann distribution of a particle with 
coordinate 770 moving in the potential field f e (r], t\t] ) 
at a temperature r. In the limit where the variation of 
the effective potential f e becomes large as compared to 
the temperature, the order parameter becomes frozen at 
the global minimum of f £ . The appearance of a global 
minimum different from 770 = is associated with a phase 
transition. Depending on whether the minimum under- 
goes a finite jump from zero or gradually moves away 
from it as a function of temperature, the phase transi- 
tion is usually said to be first or second order. 

Note that the shape of the free energy f £ as a function 
of rjo is given by a universal function independent of both 
the energy scale a and the parameter s. However, the 
amplitude of f s is inversely proportional to e. 

Free energy. Let us analyze the evolution of the free 
energy vs. 770 curve with changing temperature. Figure 
[5] corresponds to the experimental conditions of Fig. ^ 
For this set of parameters we get e = 1.77 x 1CP 4 and 
77 = 0.186. For high temperatures the only minimum of 
free energy is located at 770 = 0. At t = .0300 a sec- 
ond minimum appears and at t = t c = .0293 it becomes 
global (770 = 0.0820, A = 462), moving then gradually 
toward 7/0 = t] as the temperature decreases. Since the 
height of the free energy curve is finite, we do not expect 
any sharp change at r = r c (fceT = 1-65 hio±). Nev- 
ertheless, Fig. ^a) shows that at this temperature the 
condensate occupation does begin to increase rapidly. 

Thermodynamic limit and phase transition. We are 
going to identify a limit where the above-described quasi- 
discontinuity in the number of condensed particles be- 
comes a real discontinuity. Consider the following limit: 



FIG. 3: Phase diagram of the solitonic Bose condensate 



Notice that in this limit the Landau free energy 
/ E (?7, t I rjo) grows to infinity while preserving its shape 
along with positions of minima. A typical condensate 
occupation vs. temperature dependence in the thermo- 
dynamic limit is shown in Fig. ^b) ^^|. Figure |5] shows 
our system's phase diagram. It is universal, using the 
rescaled quantities r and 77. Our main result is the fol- 
lowing implicit equation for the critical temperature T c : 
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mgj D k B T c ) 
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Here w = 3 + W (-3e- J ) = 2.82144. . ., where W (-) 
is the principal branch of the product log (Lambert W) 
function ^ij, the inverse of x i— > xe x for x ^ — 1. Eq. (0J 
is valid in the currently experimentally accessible range 
t, 77 < 1. There r c satisfies 77 = r c ln( {3 _ 2 * )v)i ±) + 0(t?) . 
For the experiment 0, we have T c ~ 1/iK |igj|. 

Metastability. Let p No (t) be the probability distri- 
bution over A , U L (t) = YaST' 1 PM and E r(*) = 



Eil n™* P* (*) ' where £ N o 
cal maximum and nonzero minimum of f £ . For tempera- 
tures around T c and small e, the equilibrium probability 
distribution consists of a peak at Nq — (n^ 1 ) and a 
peak at A™' 11 (n^ q ). The latter dominates below T c ; the 
former, above it. In terms of the analogy to a particle 
in a potential, switching Ao from one peak to the other 
involves the system having to overcome the potential bar- 
rier at A™ ax . We assume that the two peaks satisfy the 
rate equation IIr = — IIl = j- [n^ q Hl — n^] , whose 
solution is n R , L (t) = n^ q L (l - e-*/x ) + n R , L (0) e-*/* 
Below T c , x/^r ~ X is t ne condensation time; above 
it, x/n^ q w x is the evaporation time. To obtain x-, we 
take that p N(1 satisfies a rate equation of the form p N(j = 



max, mm 



max, mm 



are the lo- 



e — > 0, r = const, 77 = const . 
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W(N *-N + l)p NQ + 1 + W(No-1^N )Pn -i ~ [w(N -l<-N ) + 



W(n —>n +i) J p NQ and that the detailed balance property 

w(N a -i^N Q )p c ^ _ 1 = w(N -l<-N )p e x holds. For sharp 

peaks we can relate the microscopic and macroscopic 

quantities through IlL/n L q ~ Pn /t>n q and Hr/ITS 5 
/ cq ..... r^—i 



This eventually yields [15| 
1 Sl Sr. 



X 



Te 2 S L + S R 



Sc , 



(5) 
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e -±fe(V,T\ V0 ) dr] r 
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and ^, T (r?o) = r | ry ). 
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Experience with Kramers' 
escape rate suggests that these integrals should be evalu- 
ated (numerically) exactly [lfij , particularly in situations 
such as that of Fig. da) when the peaks of p^ are well- 
defined but not extremely sharp. The rate coefficient 



r = w(ns 



(6) 



is the only unknown quantity, thermalization-scheme- 
dependent and beyond the scope of this paper to com- 
pute. Assuming the temperature dependence of T is 
slight next to that of x, studying Eq. shows that 
x(t) has a peak, located near r c and approaching it as 
e decreases. For parameters as in Fig. Qfa), x( T c) ~ 
70, 000 r -1 . For these parameters xi T ) does not have a 
very sharp peak, but it acquires one with decreasing e. 

Conclusions. We have presented an exactly solvable 
semi-phenomenological theory for the thermodynamics 
of harmonically trapped attractive ID bosons. Our prin- 
cipal result is that this system has a discontinuous no- 
soliton-soliton phase transition. We computed the criti- 
cal temperature and produced a phase diagram universal 
with respect to experimental details such as trapping fre- 
quencies, atomic mass, and the 3D scattering length. We 
also described the major scaling laws governing the con- 
densation time, showing it to be exponentially long as 
one approaches the thermodynamic limit. As a function 
of temperature, the transition time near the critical tem- 
perature is given by a peaked curve, which we determined 
up to a thermalization-scheme-dependent prefactor. 

Experiments thus far 0, Q make solitons by first con- 
densing repulsive atoms and then switching the sign of 
the coupling constant using the Feshbach resonance. This 
produces certain excited solitonic states that are not 
states of thermal equilibrium. In contrast, to observe 
the phase transition considered in this paper one would 
need to produce clouds of attractive ID atoms in states 
of thermal equilibrium at various temperatures. While 
challenging, experiments like this are probably already 
feasible. Conceptually the simplest strategy may be to 
perform experiments as in 0-OI- but wait long enough for 
the system to equilibrate before taking measurements. 
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One may improve this approximation by using the soliton 
energy obtained via numerical solving a Gross-Pitaevskii 
equation in presence of the harmonic trap. 
We chose parameters for this plot so that T c is slightly 
higher than the transverse level spacing. At least at the 
Boltzmann equation level, ID gas cannot thermalize by 
itself. An additional reservoir is required, thermal trans- 
verse excitations being a good candidate. In this case the 
numbers of atoms involved should be understood as those 
in the ground transverse mode of the trap. 
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[19] This temperature is in fact many times higher than the perature was much less than T c . 

transverse confinement level spacing huj±/kB, making our 
ID model invalid at T c . In reality the experimental tem- 



